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Abstract 

In this paper, we investigate the effect of adding a stabilizer on the dynamic thermal properties of ZnO nanofluid 
(containing 5 to 10 nm diameter of ZnO nanocrystals) measured using a 3oj method. Addition of the stabilizer 
leads to the stabilization of the nanofluid and also substantial reduction of the enhancement of thermal transport 
compared to that seen in the bare ZnO nanofluid. This also alters the frequency dependence of the thermal 
transport and the characteristic time scale associated with it. It is suggested that the addition of the stabilizer 
inhibits the thermodiffusion-assisted local aggregation thus leading to substantial reduction of the enhancement of 
thermal transport properties of the bare nanofluid as proposed in some recent models, and this also alters the 
characteristic time scales by altering the scale of aggregation. 
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Background 

Nanofluids are dispersions of nanoparticles (typically sizes 
approximately 5 to 20 nm) in liquid medium. In recent 
years, they have attracted considerable attention due to 
enhanced heat transport properties as seen through en- 
hanced thermal conductance [1,2]. In general, heat trans- 
port due to conducting metallic or solid inclusions in 
nonconducting fluids leads to an enhancement. However, 
in the nanofluids, which have solid inclusions of sizes in 
the range of few nanometers or few tens of nanometers, 
the enhancement in thermal conductivity was found to be 
much larger than that expected from Maxwells effective 
medium theories [3,4], 

A number of mechanisms have been proposed that 
could be responsible for the enhancement of the thermal 
conductivity. They include the (a) Brownian motion of the 
nanoparticles [5,6], (b) molecular-level layering of the 
liquid at the liquid-particle interface [7], (c) ballistic heat 
transport in nanoparticles [8], and (d) local clustering of 
nanoparticles [9,10]. The suggested mechanisms do 
provide some level of explanation of the enhancement. 
However, there is no accepted theory/mechanism that can 
explain all the observations adequately. 
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Recently reported experimental studies suggest that the 
formation of local nanoparticle aggregate can play a sig- 
nificant role in the thermal transport in nanofluids [9,10]. 
In the context of nanofluids containing Fe nanoparticles, 
it was demonstrated [11] that Fe nanoparticles in the 
nanofluids can locally assemble into aggregate of micron- 
size clusters. It was found in CuO nanofluids that large 
thermal conductivity enhancements are often accompan- 
ied by sharp viscosity that increases at low nanoparticle 
volume fractions, which has been inferred as an indicative 
of local aggregation effects [12]. The aggregation can be 
controlled by surface charge, and the critical importance 
of particle surface charge in nanofluid thermal conductiv- 
ity has been demonstrated [13]. 

In this paper, we carry out an investigation on the effect 
of local aggregation on the thermal transport in nanofluids. 
This was done in nanofluids containing ZnO nanoparticles 
with and without stabilizer. The stabilizer can affect local 
aggregation which in turn can substantially change the en- 
hancement of the thermal conduction in nanofluids. 
Importantly, we also show that this affects the characteristic 
frequency scales associated with the dynamical heat trans- 
port in such nanofluids. The characteristic frequency, as 
derived directly from the dynamic thermal property meas- 
urement, is a measure of the scale of local aggregation. The 
reduction of the scale of local aggregation can reduce the 
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magnitude of the thermal transport enhancement, provid- 
ing a direct link between the two. 

The choice of ZnO nanofluid for the investigation origi- 
nates from the fact that unlike many metallic nanofluids, 
ZnO nanofluids can be a stable suspension over hours 
even without added stabilizers. This stability arises due to 
surface charges on as-prepared ZnO nanoparticles [14]. 
The stability over hours is long enough that it enables us 
to carry out the thermal measurements. The addition of 
polyvinylpyrrolidone (PVP) as a stabilizer enhances the 
stability even further to weeks and even months. Thus, 
the system chosen is a very suitable system where the 
measurements can be carried out in nanofluids with and 
without stabilizers and thus track the changes in thermal 
parameters in the addition of the stabilizer. 

In our earlier work on ZnO nanofluids [15], which is 
carried out using a dynamic 3cj technique, it has shown 
that the parameter effusivity (Cyc, C^ = heat capacity, 
k = thermal conductivity) has a prominent frequency de- 
pendence. The measured effusivity shows appreciable en- 
hancement at low frequency, but above a characteristic 
frequency, the enhancement is significantly reduced and it 
approaches the parameters of the base liquid. In this 
paper, we investigate what happens to the enhancement of 
C p k as well as its frequency dependence when a stabilizer 
is added to the system. We find that the presence of 
stabilizer, which reduces the local aggregation, actually 
leads to a significant decrease of the C p k. We also find that 
the frequency dependence of C p i< in bare ZnO nanofluid 



gets quantitatively modified when the stabilizer is at- 
tached. In addition, we carry out an analysis of the fre- 
quency dependence of the temperature oscillation to 
separate out the contributions of C p and k components 
and find that the enhancement in C p i< is primarily due to 
the enhancement of thermal conductivity k. 

Methods 

Nanofluid synthesis 

Stable ZnO nanofluid, which is a dispersion of ZnO 
nanocrystals in ethanol, is prepared by wet chemical 
method [16]. The nanocrystals of ZnO were synthesized 
at low temperature (<90°C) in an alkaline medium using 
Zn acetate. The nanocrystals of ZnO are crystalline with 
an average size of approximately 10 nm as seen using 
the transmission electron microscope (TEM). The typ- 
ical TEM image of a nanoparticle is shown in Figure la. 
Two nanofluids were prepared. One is a pure dispersion 
of the ZnO nanocrystals in ethanol, and the other is 
made by adding PVP as a stabilizer. PVP binds to the 
polar surface of ZnO. The ZnO nanofluid, even without 
PVP, can be stabilized in scales of hours. The addition of 
PVP leads to substantial enhancement of the stability of 
the nanofluid. PVP has been used in the past to make 
stable metal colloids of Pd [17], Au, and Ag [18]. We 
have monitored the stability of the nanofluids by meas- 
uring the UV-vis spectra of the fluids as a function of 
time. Both the nanofluids show characteristic absorption 
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Figure 1 TEM image of ZnO nanocrystal used, (a) Time dependence of UV-vis spectra for ZnO nanofluids, (b) without and (c) with 
PVP stabilizer. 
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Figure 2 Typical temperature oscillation 6T 2(0 as a function of 
frequency measured in PVP-stabilized ZnO nanofluid. 



around A ~ 360 nm, which is the absorption edge for 
ZnO. For the ZnO nanofluid without PVP, the absorp- 
tion initially decreases with time as shown in Figure lb. 
The decrease is rapid initially and then slows down con- 
siderably after 30 min, when the absorption decreases by 
about 3% in 1 h. This stability is long enough to carry 
out the thermal measurements over a period of 2 h. The 
addition of the PVP leads to a very stable nanofluid that 
is stable over few weeks as can be seen in Figure lc 
where there is no perceptible change in the UV-visible 
absorption even after 2 weeks. 

Thermal measurements using 3u technique 

The thermal measurements were done using a 3co tech- 
nique [19-21], where we use a platinum film both as a 
thermometer and a heater. The method, as applied to 
nanofluids, is explained elsewhere [15]. Here, we provide a 
small gist for quick reference. 
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Figure 3 Frequency dependence of effusivity of base liquid 
ethanol, bare ZnO nanofluid, and PVP-stabilized ZnO nanofluid. 



In this method, the Pt film (width of 300 (im, thickness 
of 50 nm, and length of 5 mm grown on a glass substrate 
by magnetron sputtering) carrying a current at frequency 
/is immersed in the liquid in which measurements have to 
be made [19]. The periodic heating of the film, due to the 
sinusoidal current, makes the temperature oscillate around 
the average with an amplitude ST 2co at a frequency 2o) 
(o) = 2nf). This leads to resistance oscillations of ampli- 
tude ST 2oj at frequency 2o) around the mean, where 
SR 2co = aR 0 ST 2(O , a is the temperature coefficient of resist- 
ance (TCR) of the heater, and R 0 is the average resistance 
of the heater. The resistance oscillation SR 2oj at frequency 
2co mixes with the current at frequency o to produce a 
potential drop (V^ w ) with a component at 3co (sum band). 
The experiment measures the complex voltage V 3o) with 
its phase and amplitude, using a phase-sensitive detection 
technique. The thermal properties of the heater-on-sub- 
strate (S) and surrounding liquid (L) are given by two pa- 
rameters Z and the phase 0. These parameters are 
obtained experimentally from the observed 3co signal , 
the area of the heater (A), the power dissipated (P), and 
the measured TCR (a) of the Pt film using the equation [19] 



Zef* 



2A 
aP 



2a)V* 



(1) 



where the thermal parameter is the effusivity given as 
t\ = C p K. L and S refer to the liquid and the substrate, 
respectively. 

The Pt film has a resistance of ~ 100 £1 and a measured 
temperature coefficient of resistivity a~3.5x 10" 3 /K. The 
relative size of the heater width and the thermodiffusion 
length A t h = \/D/2co (D = thermal diffusivity) determines 
the low-frequency range of the experiment. In our case 
for the base liquid ethanol (D~9x 10~ 8 m 2 /s), the work- 
ing frequency is 1 Hz < ^ < 1 kHz for the width of the 
heater used (approximately 300 |im). At high-frequency 
range, the limit arises due to the low value of the signal. 
The experiment was carried out in a temperature- 
controlled bath stabilized to better than ±0.01 K which 
houses a cylindrical copper shell as the sample container. 
The typical data-taking time for a given frequency scan 
over the full range is 30 min. After each scan, the suspen- 
sion is shaken in an ultrasonic shaker before the next run 
begins. 

Using relation \Z s \=-^= and |Z NF | = -?= 1 7= , we 

y/$s v^ + v^NF 

obtain the ^ NF f° r the nanofluid given as [19] 



(C p k) 



NF 



t>NF 



(2) 



In addition to the effusivity we also find the ther- 
mal conductivity k using the frequency dependence of 
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Table 1 Comparison of thermal parameters for nanofluids as measured by two methods 



Quantity/method 


Bare ZnO nanofluid 


ZnO nanofluid with PVP 


Relative enhancement of effusivity £ = C p k wrt ethanol/from 3u method using V 3a) 


4.0 


2.7 


Relative enhancement of thermal conductivity k wrt ethanol/from temperature oscillation 6T 2a) 


4.2 


2.4 



the temperature oscillation ST 2oj . The ST 2oj for a line 
heater has a total width of 2b dissipating power PJxrnit 
length and immersed in a liquid [20] : 



6T : 



2co 



Pl 

TTK 



sin 2 (Kb) 



1 



{Kb) 2 Yl + Y s + YsY /interface 



-dK, 



where K is the integration variable, q 



(3) 



\J( 2 " pC >/> ' 
Yj = K jyjK 2 + qj j = s,l refer to the solid (substrate- 
carrying heater) and the liquid, respectively. The value 
of the interfacial resistance is expressed as interface- 
6.1 x 10~ 7 m 2 K/W [20]. From Equation 4, it can be 
shown that the frequency dependence of ST 2co has a 
logarithmic dependence on /whose slope is given as [21] 



2n J \d\no) J 



(4) 



We also determine the specific heat C p of the base liquid 
and the nanofluids using a differential scanning calorimeter, 
operating in modulation mode (with frequency <10 mHz). 

Results and discussions 

Change in thermal effusivity in the addition of stabilizer 

The representative data on the detected temperature os- 
cillation ST 2oj as a function of frequency is shown in 
Figure 2. It shows the typical ST 2(0 data for ZnO-PVP 



1.0 

E 

I 0.6 
O 

5 0-4 H 

Q. 

Si o.2H 

0.0 



O ZnO PVP Nanofluid 

\ 5th Order filter fit 

□ ZnO Nanofluid 
» 1st Order filter fit 




10 



f(Hz) 



100 



Figure 4 Low-pass filter response fit for ZnO nanofluids and 
ZnO-PVP nanofluid. 



nanofluids. From this data, we do the analysis of thermal 
conductivity of respective nanofluids. 

In Figure 3, we show the effusivity ^ NF = C p i< of the base 
fluid ethanol along with two nanofluids: the bare ZnO 
nanofluid as well as the ZnO nanofluid with stabilizer PVP. 
The data for the base liquid ethanol are also shown. The 
parameters are obtained from Equations 2 and 4 using the 
measured data. Both the nanofluids have the same volume 
fraction of 1.5% and have similar average particle size. 

The enhancement of ^ NF in the nanofluids, at low fre- 
quency, compared to that in ethanol is clearly seen. Im- 
portantly, it is observed that the enhancement in the 
bare nanofluid (without stabilizer) is much larger com- 
pared with that in the nanofluid with the PVP stabilizer. 
The results are summarized in Table 1, where we show 
the enhancement of the effusivity c) = C p i< as a ratio taken 
with respect to (wrt) the base fluid as determined from 
the analysis of the V 3co signal. The low-frequency-limit- 
ing values for £, were used for the parameters in Table 1. 

The value of the measured specific heat C p of the base 
fluid as well as the nanofluids are comparable (C p ~ 2.5 
J/g K). It is thus clear that the enhancement of the 
effusivity in both the nanofluids is arising primarily due 
to the enhancement of the thermal conductivity k. 

To make an independent check on the enhancement of 
the thermal conductivity, we used the measured frequency 
dependence of the thermal oscillation 8T 2o) . Equation 4 
gives a limiting low-temperature slope for ST 2(0 wrt the 
frequency (log/) that is proportional to /c" 1 . Using this in- 
formation, we obtain the relative enhancement of the ther- 
mal conductivity wrt the base fluid ethanol. The data for 
both the nanofluids are shown in Table 1. It can be seen 
that this also gives nearly the same value for enhancement 
(within 15% to 20%), which confirms that there is indeed 
an enhancement in k in the nanofluids. It is gratifying that 
the analysis from both the parameters ST 2co and gives 
similar results. 

It can be seen from Table 1 that the enhancement k 
for the bare ZnO nanofluid is significantly larger than 
that seen in the PVP-stabilized ZnO nanofluid. This 



Table 2 Corner frequency, relaxation time, and estimated 
length scale of local agglomeration obtained from the data 



Nanofluid system 

ZnO 

ZnO + PVP 



fc (Hz) 

23 ± 1 .5 
43 ± 2.3 



r(ms) 

4±3 
2±1 



La (Mm) 

18 + 2 
13 + 2 
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gives us the first important result that there is indeed a 
significant reduction in the effusivity and thermal con- 
ductivity on stabilizing the ZnO nanofluid with stabilizer 
that inhibits the local aggregation significantly, which in 
turn leads to its long-term stability. This observation es- 
tablishes a direct connection between the enhancement 
of k and the local aggregate formation. 

The frequency dependence of the enhancement and its 
analysis 

The enhancement of the effusivity in nanofluids has a 
frequency dependence as shown in Figure 3, where the 
enhancement decreased at higher frequency, and for f> 
30 Hz, the values of C p K for both the nanofluids ap- 
proach that of the base fluid ethanol. This frequency de- 
pendence of the effusivity for bare ZnO nanofluid 
(without PVP) has been reported elsewhere [15]. It was 
proposed that the frequency dependence can arise from 
dynamic local aggregation. In this paper, we explore the 
proposed hypothesis whether the frequency dependence 
indeed has a connection to the local aggregation. At low 
frequency (f< 10 Hz), the enhancement is large, and it 
reaches a frequency-independent value. 

The decrease in the effusivity at higher frequency in 
both the nanofluids can be fitted by the low-pass filter 
relation: 

The corner frequency f c and the order of the filter n 
can be obtained from the fit to the data. For the ZnO 
nanofluid without PVP, the data can be fitted by the 
first-order filter function (n = 1). For fluid with PVP, we 
got a different higher order value, which is n = 5. 

In Figure 4, we show the fit of the data to Equation 5. 
The data for both the nanofluids are shown. 

The data are summarized in Table 2. From the table, it 
can be seen that there is a significant reduction in the 
corner frequency f c in the addition of the stabilizer. 

The thermally driven local aggregation, which would en- 
hance the local thermal transport and hence the value of 
the thermal conductivity, would lead to solid-like aggre- 
gated region in the nanofluids. It is proposed that the re- 
sponse of the type shown in Equation 5 is a manifestation 
of this local aggregation. The local aggregates respond to an 
oscillating temperature field ST 2oj with a characteristic ther- 
mal relaxation time r c . This will be related to the character- 
istic length scales of the local aggregate L A through the 
thermal diffusivity D by the relation r c ~ D 1 L A . The relax- 
ation time will determine the corner frequency f c ~ {^ttt^~ 1 
(the extra factor of 2 arises because the temperature oscilla- 
tion is at frequency 2/). For frequencies larger than 2f a the 



temperature oscillation is too fast for the aggregate to re- 
spond leading to a decrease in the enhancement of heat 
transport. 

In Table 2, we show the characteristic time r c as well 
as the aggregation length L A as derived from the data. 
We find that the addition of the stabilizer leads to the 
reduction of the aggregation length L A by 25% to 30%. 
The corresponding reduction in effusivity or the thermal 
conductivity is around 40%. This agrees well with the 
hypothesis that the local aggregation can control the en- 
hancement of the thermal transport as well as the fre- 
quency response. 

Conclusions 

We have investigated the dynamical thermal property 
(effusivity and thermal conductivity) of ZnO nanofluids 
containing ZnO nanocrystals with an average diameter of 
10 nm with and without PVP stabilizer. This was done to 
investigate the role of the stabilizer in the enhancement of 
thermal transport properties of nanofluids. It had been 
suggested that thermodiffusion-assisted solid-like' local 
aggregation of the nanoparticles in the nanofluids can be 
the origin of enhancement of thermal conductivity in 
nanofluids. The investigations carried out on bare ZnO 
nanofluids as well as PVP-stabilized nanofluids show that 
addition of a stabilizer, which inhibits diffusion-assisted 
local aggregation due to attached moiety, leads to reduc- 
tion in the enhancement of thermal parameters that are 
observed in bare ZnO nanofluids. It has also been shown, 
from characteristic time scales of the dynamic thermal 
measurements, that the scale of aggregation gets reduced 
in the addition of stabilizers. The experimental results 
provide evidence that the origin of enhancement of ther- 
mal conductivity in nanofluids can arise from local aggre- 
gation that occurs by thermodiffusion. 
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